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Abstract — A recent paper [1] (El Aroudi, 2012) misapplied 
a critical condition (Fang and Abed, 2001) to a well-known 
example. Even if the mistake is corrected, the results in [1] 
are applicable only to buck converters and period-doubling 
bifurcation. Actually, these results are known in Fang's works a 
decade ago which have broader critical conditions applicable to 
other converters and bifurcations. The flaws in [1] are identified. 

Index Terms — DC-DC conversion, subharmonic oscillation 



I. Introduction 

Recently the author of [111, based on ||2l, presents critical 
conditions for current/voltage mode control (CMC/VMC) for 
buck converters. The applicability of critical conditions is not 
stated, which may mislead or confuse the readers. As an author 
of [2], I feel a need to clarify the results in [1]. Actually, the 
results are are known in 0, H, ||5], dU, d, li], S a decade 
ago which have broader critical conditions applicable to other 
converters and bifurcations. The author of [1] was aware of 
these results when he received an email copy of [3] on Aug. 6, 
2010 with the key results particularly identified, but [3] was 
not cited in [1]. 

II. Broader Results Known a Decade Ago 
A. Unified VMC/CMC Model With a Common Ramp 

In VMC, a sawtooth signal is used for PWM modulation. 
In CMC, a ramp signal is used for stabilization. The two 
signals serve different purposes. However, they have the same 
waveform. In ||3|, ID, ||5|, ||6|, they are modeled by the same 
signal h{t) in a unified VMC/CMC model shown in Fig. [T] 
where h{t) := Vi + {Vh — Vi){-^ mod T) varies from a low 
value V; to a high value Vh with an amplitude of Vm = Vh — Vi. 

Since the model is applicable to both VMC and CMC, the 
system dynamics and critical (boundary) conditions associated 
with the unified model are also applicable to both VMC and 
CMC. No special or separate analysis is required, and the 
unified model has broad applications. 

Also, the unified model is applicable to both the power stage 
and the closed-loop converter (power stage plus compensator). 
For the power stage, a control signal is used to control the 
inductor current in CMC or the output voltage Vq in VMC. 
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Figure 1. Unified VMC/CMC model for a switching conveiter. 

For the closed-loop converter, a reference signal is used for 
similar purposes. The two signals are modeled as iv. Let the 
source voltage be Vg. Both v^ and v^ are external signals input 
into the converter, collectively represented by a vector u. 
In the model, Ai,A2 G R^'^^, Bi,B2 £ R^^^ 



C,Ei,E2 e R^^'\ and Z? G R 



1 1x2 



are constant matnces. 



where N is the system dimension. Within a clock period T, 
the dynamics is switched between two stages, 5*1 and S2. Let 
the compensator output be y := Cx + Du e R. The dynamics 
is switched to Si at t = nT. The dynamics is switched to S2 
when h{t) > y{t). In the trailing-edge modulation (TEM), Si 
is the ON stage, and 5*2 is the OFF stage. In the leading-edge 
modulation (LEM), Si is the OFF stage, and S2 is the ON 
stage. Denote the switching frequency as fg := 1/T and let 
Ws := 271 fs. 

B. Closed-Form Jacobian Matrix, First Known in 

The periodic solution x^{t) of the system in Fig. [T| corre- 
sponds to a fixed point x^{0) in the sampled-data dynamics. 
Let x°{d-) = Aix°{d) + Biu and x"{d+) ^ A2X°{d) + B2U 
denote the time derivative of x'^{t) at t = d^ and 
respectively. Let y'^{t) = Ca;°(t) + Du. In steady state, 
= Ci°(t). Let the steady-state duty cycle be D and 
let d be the switching instant within a cycle. For LEM, 
d^{l- D)T; for TEM, d = DT. Confusion of notations for 
capacitance C and duty cycle D with the matrices C and D 
can be avoided from the context. 
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In steady state, 



(1) 



c°(0) 



T-d 



= Ca;"(d)+Du = /i(d) 



e^^''daB2U (2) 



(3) 



Solving ([TJ-(l3]i, one obtains and x^{d). 

Using a hat " to denote small perturbations (e.g., Xn = 
Xn — a;°(0)), where x„ is the sampled state at t = nT. From 
El, H, Is), the linearized sampled-data dynamics is 

Xn+l = (4) 

where the Jacobian matrix, ^rif known in 13], ID, is 

CiO(d-)-/i(d) 

A closed form like Q is important, because it leads to many 
general critical conditions discussed next. 

C. Equivalent Linearized Sampled-Data Dynamics 
Let $ = $0 - rvp, where 

$0 = ^A2{T-d)^A^d 



r 



eM^-''\x''{d-)~x\d+)) (7) 



(8) 



Cx°{d')-h{d) 

The dynamics (|4|l can be transformed into a plant (having 
an input dn and an output ?i„ = with a Mn/fy negative 

feedback {dn = —Wn), 



Xn+l = ^oXn + rd,, 
dn = -'^Xn 



(9) 



D. Broader General Closed-Form Critical Conditions 

Theorem 1: ||3] p. 46] Suppose that A is not an eigenvalue 
of <i>o. Then A is an eigenvalue of $ if and only if 



Cx^id-) + Ce'^i'*(A/ - $o)"'r = hid) 



(10) 



Proof: Suppose A is not an eigenvalue of $o. 

dot [A/ - = dct [A/ - $o] dot [/ + ( A/ - $0 ) " ^ r*] 
= dct[A/-$o](l + *(A/-$o)"^r) 

Then dct[A/ - $] = requires that *(A/ - $o)"^r = -1, 
which leads to ( fTOb by using (|6]l-(l8]l. □ 
Instability occurs when there exists an eigenvalue A (also the 
sampled-data pole) of $ outside the unit circle of the complex 
plane, which leads to the following corollary, associated with 
three typical instabilities [3], Q, E): period-doubling bifurca- 
tion (PDB, subharmonic oscillation), saddle-node bifurcation 
(SNB), and Neimark-Sacker bifurcation (NSB). 



Corollary 1: 13] p. 46] (i) If the system parameters corre- 
spond to an occurrence of PDB (A = —1), then 

Cx'^{d-)-Ce^"^{I + <^>o)-^r^h{d) (11) 

(ii) If the system parameters correspond to an occurrence of 
SNB (A = 1), then 

Ci°(d-) + Ce^i'^(/-$o)"^r = /i(d) (12) 

(iii) If the system parameters correspond to an occurrence of 
NSB (A = eJ", 61 7^ or tt), then 

Cx°{d-) + Ce^^'^ie'^I - <Po)-^T = h{d) (13) 

All of the results above are applicable to general converters, 
they are applied to a special case: buck converters. 



E. Applying General PDB Condition i liiD to Buck Converter 

Let Bi := [Bii,Bi2], B2 := [B2i,B22] to expand the 
matrices into two columns. Many LEM and TEM examples 
are presented in ||3]. 

The LEM buck converter generally has 

Ai=A2=A, Bi = [0Arxl,Bl2], B2 = [B,Bi2] (14) 

Using ([T]i and (|2]l, the PDB critical condition ( fTTT i becomes 

C[{I+e-^^)-' + {I-e^^)-\e''^ -e^')]Bv, = h{d) (15) 
or in terms of Vs, shown in 15] Eq. 3.118], with d ~ {1 — D)T, 

h{d) 



C[{I + e-AT)-l + (J _ gAT)-l(gAT „ e^T(l-Z5)]5 

(16) 



F. Directly Applying (|77} to TEM: Different Dynamics 
The TEM buck converter generally has 

Al=A2-A, Bi = [B,Bi2], B2 - [0Arxl,Bl2] (17) 

Using ([T]i and (|2]l, the PDB critical condition ( fTTT i becomes 

C[{I - e^^)-i(e^'^ + + e^^)-^]Bvs = h{d) (18) 
or in terms of Vs, with d = DT, 

" ^ ^ ~ C[(/-e^^)-i(e^^^- /) + (/ + 6-4^)-!]^ 

(19) 

It is known decades ago that TEM and LEM have different 
dynamics Qol, CHI- Thus, vJ^^{D) 7^ v^^^{D). However, 
[1] wrongly assumes vJ^^{D) = v^^^{D). As a side note, 
the recent results in lfT2l . lfT3l were repeatedly said to be wrong 
by a reviewer on this issue. That reviewer was wrong. One can 
prove vJ^^{D) = -v^^^{l - D) by simple algebra. 
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G. Equivalent Critical Condition Based on Harmonic Balance 

Let Re denote taking the real part of a complex number 
Let the voltage across the second switch (or the diode) be Vd- 
Let the Vd-to-y transfer function be G{s) for LEM, or -G{s) 
for TEM. Based on harmonic balance analysis, an equivalent 
PDB critical condition for the buck converter (due to its special 
property ([14) or ([17}) is obtained |l2], ||3], JHl, 



2Re 



.fc=i 



(20) 

As reported in ||3] pp. 72-73], both (|20) and ([16) are exact 
PDB conditions, and they are the same. Therefore, 



2/.Re 



.k=l 
AT\ 



C[{I + e-^'^')-i + (/ - e^^)-\e^^ - e^'^)]B (21) 

Note that dTSl l. dTSl l. (i20] i. and dlTT i. expressed in terms of 
d, are valid for both LEM and TEM (d = (1 - D)r for 
LEM, d = DT for TEM), whereas the main result of [1], (17) 
(underline added for [1]), is limited to TEM only. 



III. Simple Direct Extension of (fTSl) 

Using Taylor series expansion, ([TSl l directly leads to, also 
shown in |[T2l submitted seven months earlier than [1], 



C{Y^5^{D)A^T^)Bv, ^ h{d) 



(22) 



n=0 



where 6o{D) = (1 - 2L»)/2, (5i(i:>) = (-1 + 2D- 2L>2)/4, 
and 52{D) = {-D + - 2D^)/12. If the value of a real 
pole (eigenvalue of A) is comparable with ujs, the high order 
terms in ( |22] | are significant and cannot be ignored |[T2l . 
Note: (17) of [1] is exactly ([lUl, and (20) is exactly (l22T i. 

IV. Flaws, Mistakes, and Confusions in |[T1 

First, all of the obtained results in [Tj are applicable to 
buck converters only, and therefore the title of |[T] without this 
limitation is misleading. 

Second, the condition ( i20l ) is applied without stating its 
applicabiUty. In d = (1 - D)T for LEM. However, 

d = DT is used instead in [1]. The mistake occurs at 
the beginning and all the results in [JJ depend on proper 
application of ( |20] |. 

Third, throughout [T], different TEM and LEM examples 
are commingled. It is probably assumed in yj that any critical 
condition is applicable to both, which is wrong. If the goal of 
[1] is on the TEM conditions, then a wrong TEM condition is 
misapplied to the well-known LEM example on page 5. If the 
goal is on the LEM conditions, in contrary, all results (5)-(24) 
of [1] are valid only for TEM. The sequence of presentation 
in [1] is confusing: an LEM condition is misapplied with d = 
DT (actually for TEM), then some TEM conditions are derived 
but applied to an LEM example on page 5. 



Fourth, readers are misled about the cited references. On 
page 2, it is stated, "in [14], the series S{D) in (5) has 
been approximated by the term that involves the transfer 
function H{s) with the smallest argument or graphically 
solved by truncating the series to an arbitrarily large number 
of harmonics. In this work, a closed form expression will be 
given for the series S{D) defined in (5)." It sounds Uke a 
different approach is used in [1]. Actually, still a series closed 
form (20) is presented, and eventually still a figure (Fig. 3) is 
used to obtain the critical value. In [14], the approximation is 
just one way to obtain the critical value. The key in [14] is 
on the exact condition (5), on which all of the results of |[T] 
are based. 

Fifth, great efforts (from (5) to (15)) are made in [1] to 
paraphrase the equality i2T[ already reported in f3} pp. 72-73]. 
As discussed above, (15) for the buck converter is just a special 
case of the known condition ( ITOl l. Note that the condition ( ITOb 
is not only broader but its derivation is also simpler, without 
the need of the great efforts in [[jj. 

Sixth, the author of [1] is aware of the general condition 
dlOl l available in |(3] and the unified WMC/CMC model E), id, 
Ha, 0, Q, i] of Fig. [U However, none of 0, gl, E, W\, 
||7l, ID is cited. Before 1997, to my knowledge, the sawtooth 
signal in VMC and the stabilization ramp in CMC were not 
modeled by a single signal. 

Seventh, some parts are inconsistent with the others in [1]. 
Using (21) does not lead to Table I in [1], unless (21) is scaled 
by l/(27r)™. Even with this corrected, using Table I does not 
lead to (4), unless Si{D) = D — 1/2 is corrected in Table I. 

Eighth, the TEM condition (17) is misapplied to the well- 
known LEM example on page 5. 

Ninth, the readers are not able to reproduce the same results. 
In Fig. 3 of m, a plot is falsely claimed to be based on (17), 
but actually based on dTSl l already known in [3]. As stated 
in Q, gC = (8.4, 0) for this example. The plot of (17) is 
actually all negative (not shown). If it is corrected with gC — 
(—8.4, 0), the plot of (17) as shown in Fig. [2] is still incorrect 
because wrong v* ~ 35 is obtained, not v* = 31 as claimed 
in Hj. As expected, (17) is for TEM, not applicable to this 
LEM example. The correct ([TSl l for LEM (known in |[3l but 
never obtained in [1]) shown in Fig. [2] produces the correct 
V* = 31. As expected, due to vJ^^{D) = -vY''^{\ - D), 
the correct plot is symmetric to the wrong one with respect to 
D = 1/2. The plot in Fig. 3 is not based on (17) as claimed. 
How the plot was produced to obtain v* = 31 is a mystery. 

Tenth, no experimental prototype is actually shown. 

V. Conclusion and Recent Developments 

Compared with [1], broader results ||3],|i],||5],||6],||7],|[8], 
||9l known a decade ago are presented in this note. Actually, 
|[T] is inconsistent with the past works of its own author In the 
last six years, the author of [1] advocated a ripple index |[T5l 
to predict PDB for VMC. As indicated in [1], the author of [1] 
still up to now does not distinguish TEM from LEM [known 
with different dynamics), this gives doubt to the validity of his 
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Figure 2. Dashed line based on tlie claimed (17) which leads to wrong v*\ 
solid hne based on the correct M6\ known 15 years ago ||3] Eq. 3.118]. 

similar past works, further indicating that the ripple index has 
limited applicability. Actually using (5), which totally relies 
on, many counter-examples of the ripple index hypothesis can 
be obtained lfT3l . |[T6l . In fact the author of [1] attended the 
same conference session a decade ago when a counter-example 
(for average current control known to resemble VMC IfTTl ) 
was presented ifTSl . The author of [1] was fully aware of 
this counter-example, never cited ifTSl . and still advocated the 
ripple index in the last six years. 

The readers are advised to refer to some direct extensions 
lEl, m, HE, ESI, EQI, IE] based on the key general 
condition ( fTOl l. Let the left side of ( fTOl i be called an "S-plot" 
S{X, D), a function of A and D, for example. It predicts the 
occurrence of three typical instabilities (PDB, SNB, NSB) in 
DC-DC converters and shows the required ramp slope for 
stabilization. Let the loop gain transfer function of (|9) be 
7V(z). One has J\f{z) = *(z/- (f>o)"iro. The (discrete-time) 
Nyquist plot is M{e^'^'^) as a function of uj swept from to 
ujs. Also, define an "F-plot" as F{9) = S{e^^,D) swept from 
— TT to TT. In a single plot in the complex plane, either Nyquist 
plot or F-plot also predicts the three typical instabilities. 

The broader general critical conditions (fTOt-lflJll have not 
been published except in my thesis, available online since 
1997. These important and subsequent results |[T2| . |[T3l . |[T6| . 
II19I . II20I . II21L much broader than [J^ and applicable to other 
converters and bifurcations, and variable switching frequency 
case were fiercely opposed after submitted by a reviewer who 
was actually wrong about LEM/TEM. Only a small subset 
of key results (resubmitted to other journals) is about to be 
published ED, ||23l. Many key results based on (fT0li-(fT3Tl. 
although available since 1997, remain unpublished. 
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